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In the developmental phase, the sensitivity of the model 
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investigated. The results of this work were used to determine 
the optimal time step and model domain size for the main 
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In the experimental phase, the model was spun-up using 
time averaged wind stresses. The model ocean was then exposed 
to two years of realistic daily varying wind stresses covering 
the period of 1991 and 1992. The model developed an El Nino 
like response that corresponded in several respects with 
observed features of the 1991-92 El Nino. The model also 
developed tropical instability waves in the eastern Pacific 
similar to those observed in situ and in satellite SST images. 
The model's responses to the tropical cyclones that occurred 
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I . INTRODUCTION 

A. EL NINO AND WESTERLY WIND BURSTS 

The El Nino/Southern Oscillation (ENSO) phenomenon has 
received much scientific interest due to its global impact on 
the oceans and the atmosphere. The Tropical Ocean Global 
Atmosphere Coupled Ocean Atmosphere Response Experiment (TOGA- 
COARE) has as one of its major goals to provide additional 
detailed information to better understand the air-sea 
interaction in the tropics and ENSO (Webster and Lukas 1992). 
Although much is known about ENSO events and the changes that 
take place in the Pacific during these events, no single 
theory can yet explain all of the factors leading to an El 
Nino or what causes the ocean to return to its pre-El Nino 
state . 

Prior to the onset of an El Nino event, the atmospheric 
sea level pressure gradient weakens between Darwin and Tahiti, 
with the pressure at Darwin rising and the pressure at Tahiti 
falling. This pressure change is measured by the Southern 
Oscillation Index (SOI). El Nino events are associated with 
negative SOI values. A change in the atmospheric sea level 
pressure gradient can decrease the trade wind easterlies 
allowing warm water to flow to the east, which contributes to 
large scale warming in the eastern Pacific. However, we can 



1 



also follow the opposite path to get to the same result. Warm 
water flowing to the east can cause warming which will change 
the pressure gradient and weaken the trade wind easterlies. 
This feedback mechanism between the ocean and the atmosphere 
makes the prediction, modeling, and understanding of El Nino 
extremely difficult. Obviously, the complex interactions 
between the atmosphere and ocean are of crucial importance to 
the understanding of the problem (Philander et al. 1989). 
This study focuses on the importance of the forcing applied to 
the ocean by the wind. 

Several theories have been proposed on the role that 
westerly wind bursts play in the weakening of the trade winds 
(e.g., McPhaden 1988; Wyrtki 1975; Giese and Harrison 1991; 
Philander 1981). Keen (1982) proposed that sea surface 
temperature (SST) changes in the tropical Pacific at the start 
of an ENSO create favorable conditions for the formation of 
westerly winds and, specifically, for the development of 
tropical cyclone pairs that lead to large westerly wind 
bursts. Giese and Harrison (1990) suggested that a series of 
westerly wind bursts timed to have an additive effect could be 
the trigger for an El Nino and/or help to maintain an El Nino. 
Thus, westerly wind bursts may play an important role in both 
the onset and maintenance of El Nino events. 

Definitions of a westerly wind event or burst vary. Riehl 
(1979) used the term to describe anomalous low-level synoptic 
scale westerly winds. This definition includes the relaxation 
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of normal easterly flow. Thus, for example, a sudden decrease 
from normally strong easterly flow to weak easterly flow would 
be considered a westerly wind burst (Luther et al. 1983). I 
will use Riehl's definition in this paper. 

Numerous sources of westerly wind bursts in the tropics 
have been found. If, during the northern winter, the Siberian 
anticyclone shifts south, sea level pressure in the tropical 
western Pacific may rise and lead to westerly winds (Chu 
1988). A surge of the trade winds in the northern hemisphere 
that pushes the intertropical convergence zone (ITCZ) south of 
10°S can lead to westerly winds at the equator (Riehl 1979). 
Also, the formation of a pair of tropical cyclones, with one 
north and one south of the equator, can cause intense westerly 
winds for extended periods of time (Keen 1982). A single 
tropical cyclone can also cause westerly wind bursts and may 
cause westerly wind bursts that precede or lag the time the 
storm intensifies to strength T4 (32 m/s) by up to seven days 
(Harrison and Giese 1991). 

Westerly wind bursts occur frequently in the tropical 
Pacific. Rasmusson and Carpenter (1982) have developed a 
classification scheme of westerly wind bursts based on the 
latitude of the wind burst. Their study showed that there is 
an increase in the number of events just prior to and during 
an El Nino event and that most of the wind events at the start 
of an event are either centered on or just south of the 
equator. They found that wind events to the north of the 
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equator, or to the far south of the equator, are not as 
closely linked to the start of an El Nino. 

Westerly wind bursts are probably not sufficient in and of 
themselves to cause an El Nino. Wyrtki (1975) proposed that 
the Pacific must be preconditioned in order for an El Nino 
event to occur. He found that before most El Nino events 
there was an increase in strength of the southeasterly winds 
for up to two years. These stronger winds would increase the 
amount of warm water piling up in the tropical western Pacific 
and set the stage for the warm water to surge back to the 
east. This combination of a preconditioned ocean with strong 
westerly wind bursts is more likely to cause an El Nino than 
the westerlies alone. 

Philander (1981) suggested that the warming of the 
equatorial central Pacific during an El Nino event causes 
westerly wind bursts by weakening the easterly winds. These 
changes in the winds will tend to maintain the event. 

Westerly wind bursts may initiate or enhance an El Nino by 
several processes. In an observational study. Cooper (1992) 
found that a westerly wind burst can have a local as well as 
a remote response. Locally it can cause a reversal of the 
South Equatorial Current (SEC), downwelling, and a deepening 
of the thermocline. Some extremely strong events can even 
cause an eastward surface current (a Yoshida jet) to develop 
(McPhaden et al. 1988). The remote response is in the form of 
an equatorially trapped Kelvin wave that can travel across the 
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entire Pacific (e.g., Delcroix 1991, Cooper 1992). These 
eastward currents and the currents associated with the Kelvin 
wave can lead to the advection of warm sea surface 
temperatures (SSTs) thus contributing to the warming in the 
tropical eastern Pacific (Harrison and Schopf 1984). These 
Kelvin waves also cause a rise in the ocean surface and a 
deepening of the thermocline. 

Due to the small amount of data that has been collected in 
the tropical Pacific that can be used to study El Nino, models 
are needed to fill in the gaps. These models can not only 
provide insight into what is observed in the oceans, they can 
also point out areas where more observation is needed. 

B. NUMERICAL MODELING IN THE TROPICS 

A wide variety of models with varying complexity have been 
used to simulate the ENSO phenomenon. The earliest attempts 
to model the ENSO phenomenon used very simple conceptual 
models. These types of models are useful for determining the 
types of waves and instabilities that may play a role in El 
Nino. McCreary and Anderson (1991) provide a good review of 
these models. 

Models of intermediate complexity have also been useful. 
These models use detailed physics but are limited in some way. 
For example, they may use a limited domain, or they may 
prescribe the atmosphere in terms of specified atmospheric 
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parameters, and then evaluate the ocean's response as these 
parameters are varied. 

The most complex kind of model is the coupled global 
circulation model. These models combine an atmospheric model 
with an ocean model and use the solutions from each model to 
force the other one. Presently many problems still exist with 
coupled models. These problems include climate drifts that 
slowly cool the atmosphere and the ocean, and nonlinear 
instabilities. For a good review of several coupled models 
and the problems inherent in them, see McCreary and Anderson 
(1991) or Charnock and Philander (1989). 

To lay the groundwork for the work done in this paper, 
this review will concentrate on work with ocean models having 
specified atmospheric forcing. These models have received 
considerable attention and have greatly enhanced our 
understanding of the ENSO process. 

McCreary (1976) used a simple analytical model with two 
layers to simulate the tropical Pacific circulation. Although 
this model did not attempt to accurately duplicate the 
structure of the Pacific, it did show that changes in the 
winds near the equator caused a greater response than do winds 
some distance away (±5° of latitude). The model also was able 
to produce a response which resembled an El Nino in many ways. 
Later efforts with the same model, after some improvements, 
showed that if the winds are linked to the sea surface 
temperature so that a warming in the east weakens the winds, 



6 



the model oscillated on time scales similar to observed El 



Nino scales (McCreary and Anderson 1984). An example of 
attempts to use more realistic physics in the modeling problem 
can be seen in Semtner and Holland's ( 1980) work using a 
primitive equation model with 14 levels. Although this model 
was limited in basin size and had a constant wind stress, it 
did develop many of the features of the tropical Pacific 
ocean. 

The next step toward investigating El Nino came when these 
models were exposed to realistic forcings. Latif (1986) ran 
a primitive equation model forced with real wind data for 
several decades. Although the wind fields he used had all 
been smoothed or were interpolated from monthly mean values, 
he did obtain relatively realistic simulations of the 
climatologic currents and temperature fields for the Pacific. 
The model was also used to simulate the 1982-83 El Nino. The 
model developed an El Nino-like response, with tropical 
eastern Pacific SST anomalies as large as 4°C. A weakened 
equatorial undercurrent and eastward surface currents were 
also produced. Philander et al. (1987) used a primitive 
equation model of the tropical Pacific with a 33 km meridional 
grid near the equator increasing to 200 km at about 27° (the 
northern boundary was at 50 °N and the southern boundary was at 
28°S) and a fixed 100 km spacing in the longitudinal 
direction. This model included land masses if they projected 
into the area of interest; otherwise, solid boundaries were 
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used to close the domain. The model was integrated for 
several years with monthly mean winds. The model developed 
all of the major currents of the Pacific equatorial region and 
showed the seasonal dependence of these currents. 

Miller et al. (1992) used a primitive equation model in 
isopycnic coordinates with an imbedded bulk mixed layer model 
and prescribed climatologic monthly mean atmospheric 
temperature, humidity, cloudiness, and turbulent kinetic 
energy in the form of u* (surface friction velocity). The 
wind stresses were smoothed monthly means interpolated to 
daily values. Although the model did develop the major 
currents and was within one degree of the real sea surface 
temperature for most of the domain, the speeds of the currents 
were only about 50% of those observed. 

Whitney (1992) used an ocean model with an imbedded mixed 
layer model to evaluate the importance of detailed modeling of 
the upper layers of the ocean. She exposed her model to a 
highly smoothed time varying wind to evaluate the model's 
response to westerly wind events. Her results showed that 
detailed modeling of the mixed layer of the ocean is very 
important in correctly modeling El Nino. 

These models have succeeded in simulating several, but not 
all, of the major ocean features of El Nino. 
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C. UNRESOLVED PROBLEMS 

Most of the modeling efforts to date have looked at the 
responses to long term or monthly changes in the winds, while 
the observational studies (Cooper 1992, Chen et al. 1988) 
indicate that the daily and weekly fluctuations of the wind 
can play an important part in the structure of the ocean 
currents. 

The focus of this study is on the effects of the daily to 
weekly variations in wind forcing and the role they play in an 
El Nino event and in the development of the temperature and 
current patterns in the tropical Pacific ocean. In this 

study, I test the hypothesis that the inclusion of the short 

*» 

term variability in the winds over the Pacific is extremely 
important in the development of a realistic El Nino response, 
and in developing a realistic model of the tropics. 

Since the effect of the wind is initially at the surface, 
the model used for this study is a general circulation model 
with explicit mixed layer physics to enhance the modeling of 
the changes due to wind forcing. This model was first used by 
Whitney (1992) to examine the response of the tropical Pacific 
Ocean to fixed and variable winds. The relatively simple 
atmospheric forcing used in this model allows for the 
isolation of independent forcing terms to study their 
individual effects on the model ocean. 

The wind stresses used for this study were developed from 
the U. S. Navy's Operational Global Atmospheric Prediction 
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System (NOGAPS) daily average surface wind analyses without 
any spatial smoothing. In the final model runs, the model was 
initially spun-up with a two-year time mean wind stress 
calculated from the daily NOGAPS winds for 1 January 1991 
through 31 December 1992. The model was then run out for two 
years using the actual daily average wind stresses, in an 
attempt to model the 1991-92 El Nino. 

The goals of this study are: 

1. To further refine the model first used by Whitney (1992) 
to study the tropical Pacific. 

2. To test the model with realistic daily wind stress 
forcing. 

3. To examine the model's waves and currents created by 
westerly wind bursts. 

4. To perform sensitivity studies to identify ways to 
improve the model. 

Chapter II contains a description of the model and the 
wind stress fields used in this experiment. Chapter III 
contains the results of some of the developmental and the 
experimental model runs. Chapter IV is a discussion of the 
results of the model runs, some conclusions, and suggestions 
for further improvements to the model and areas needing 
further investigation. 
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II. 



THE MODEL 



A. MODEL FUNDAMENTALS 

The model used for this study is a 15-level primitive 
equation ocean general circulation model (OGCM) with an 
imbedded bulk mixed layer model. The model was originally 
developed in radial coordinates and used to study the ocean 
response to a hurricane (Adamec et al . 1981). The model has 
since been reformulated in cartesian coordinates and has 
proven to be extremely versatile. Whitney (1992) adapted the 
model for use in the equatorial region and performed much of 
the initial developmental work for application to the tropical 
oceans. The model for this study is basically the same as 
that developed and tested by Whitney (1992). However, a few 
minor modifications have been made to improve the dynamics and 
response of the model. 

The OGCM portion of the model was originally developed by 
Haney (1974) as a 6-level primitive equation model used to 
study oceanic responses to large-scale heat and momentum 
fluxes . 

The mixed layer model was developed by Garwood (1977) and 
is based on the turbulent kinetic energy equation. The model 
assumes that the water column below the mixed layer is 
dynamically stable, and the mixed layer is dynamically 
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unstable. Thus, the mixed layer depth is the depth to which 
turbulent mixing penetrates. The temperature and horizontal 
velocity in the mixed layer are assumed to be vertically 
homogeneous. Mixed layer conditions are allowed to change due 
to surface fluxes, advection, and entrainment. 

Adamec et al. (1981) incorporated the Garwood mixed layer 
model into Haney's OGCM to give an OGCM with realistic surface 
layer physics. Interaction between the fixed level OGCM and 
the mixed layer portion occurs in two steps. First, the 
dynamic portion of the model calculates advective and 
diffusive changes for the mixed layer model. The mixed layer 
model then computes changes in mixed layer depth, velocity, 
and temperature due to surface fluxes and entrainment and 
passes the information back to the dynamic part of the model. 
For more information on the coupling of the two models, see 
Adamec et al . (1981). 

The basic assumptions of the model are as follows: 

1. The ocean is hydrostatic and incompressible. 

2. Density is a function of temperature only. 

3. Vertical velocity is zero at the ocean surface and at the 
bottom. 

4. All of the boundaries (other than the ocean surface) are 
treated as free slip surfaces. 

5. Heat is allowed to diffuse through the north and south 
boundaries . 
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These conditions do not allow a barotropic mode. Laplacian 
diffusion is used in the horizontal momentum and temperature 
tendency equations to simulate subgrid scale horizontal fluxes 
of heat and momentum. Vertical turbulent fluxes below the 
mixed layer are also parameterized in the same manner. Table 
1 gives the values of the model coefficients. 

The model domain is rectangular in the horizontal and 
vertical, with the longitudinal extent exceeding the 
latitudinal extent in all cases for this study. The model is 
symmetric about the equator, with the middle line of grid 
points located on the equator. The 15 levels in the vertical 
are concentrated in the upper ocean to provide for better 
resolution of the mixed layer and seasonal thermocline. All 
of the boundaries are considered to be solid (no flow normal 
to the boundary) and flat (no topography) . Figure 1 provides 
a schematic of the model domain and shows the distribution of 
levels in the vertical. The upper level of the model is at 10 
meters; however, this level will be treated as representing 
the sea surface temperature (SST) . This is not a bad 
assumption, except in the case of rapid heating of the surface 
or a mixed layer depth shallower than 10 meters. 

Several different model domains were used in the 
developmental stages and in the final experiments. Figure 2 
shows the location of the different model domains in the 
Pacific and the grid sizes and number of grid points in each 
direction for each of the model runs. Five different model 
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domains were used during this study, the smallest of which 
extends from 10 °N to 10°S and covers 30° in longitude from 
170 °W to 140°W. The next model had twice as many grid points 
in both directions and extended from 20 °N to 20 °S and covered 
approximately 60° of longitude from 175 °E to 125 °W. A larger 
version of this model was used that had the same latitudinal 
extent as the previous model but had a larger longitudinal 
grid spacing so that it covered 120° of longitude from 150 °E 
to 90 °W. A later version of this same model was used with 15 
more grid points in the longitudinal direction to increase the 
basin length to 135°. This model was also centered further to 
the west than the 120° wide basin, with its western boundary 
at 135 °E and extending to 90 °W. The final model domain 
extended from 45 °N to 45 °S and covered 90° of longitude from 
160 °E to 110 °W. This domain used a variable meridional grid 
spacing to achieve the larger north south boundaries. These 
domains are referred to as the small, medium, large, large2, 
and variable grid domains in Figure 2. 

The model is initiated from rest with the temperature 
profile uniform in the horizontal and a vertical profile as 
illustrated in Figure 3. This temperature profile was chosen 
to match the profile used by Semtner and Holland (1980) so 
that the initial model runs could be compared to their work. 
The model is integrated forward in time using a leapfrog 
scheme with an Euler backward scheme applied every eleven time 
steps to prevent the solutions from separating in time. The 
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finite is differencing done on an Arakawa B grid (Haltiner et 
al. 1980) . The grid spacing, time step, and other parameters 
are given in Table 1. 

The two surface forcings required for the model are wind 
stress and the surface heat flux. The surface heat flux is 
computed using the method of Haney (1971) and is proportional 
to the air-sea temperature difference, where the net downward 
heat flux is given by: 

Q=y (T a -T s ) 



Where 

y = heat flux coefficient (W m' 2 "C' 1 ) 

T a = air temperature (°C) 

T s = sea surface temperature (°C) 

All of the model runs in this study used a fixed air 
temperature of 25 °C over the entire domain. The wind stresses 
used for the experimental model runs are described in detail 
in the next section. For the developmental runs, a constant 
westward wind stress of 0.6 dynes/ cm 2 was used. 

B. WIND STRESS FIELDS 

Wind stress is the only forcing term that was varied while 
the model was running. Initial model development was carried 
out using a spatially and temporally constant wind stress to 
provide a baseline for comparison with previous efforts (e.g., 
Semtner and Holland 1980) . This also helped in finding 
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numerical instabilities in the model (cf. Whitney 1992). 
Next, wind stresses that varied in space but not in time were 
used to produce a temperature and current pattern that more 
closely simulated the equatorial Pacific. Whitney (1992) 
added an observed time varying wind stress anomaly to a 
constant analytic wind stress field to test the model's 
response to time varying wind stress. The results she 
obtained were encouraging and led to the current work. For 
this study, we applied full-strength time varying wind 
stresses to the model. 

The wind stresses chosen for this study were derived from 
daily averages of the NOGAPS analyzed 10 meter wind fields for 
1991 and 1992. Goerss and Phoebus (1993) describe the NOGAPS 
analysis procedure. The original data set, provided by Dr. R. 
Gelaro at the Naval Research Laboratory (NRL) in Monterey, 
California, contained global wind values at 10 m above the 
surface twice daily at every 2.5 degrees of longitude and 
latitude. These winds were averaged in time to provide single 
daily average values. The winds were then linearly 
interpolated to the various model grids. The wind stress was 
calculated from these winds at each grid point using the 
formulas: 




where 




p a = the density of air (1.2 kg-m' 3 ) 
C D = drag coefficient (1.2 xlO' 3 ) 
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u 10 and v 10 = the u and v components of the wind at 10 m 
above the surface in (m/s) 

The wind stresses were then converted to dynes/ cm 2 . 

The wind stresses were also averaged over the entire two- 
year period to provide a mean wind stress field for use in 
initializing the model during spin-up. Several of the early 
model runs were performed using an average wind stress field 
created from the December 1991, January 1992, and February 
1992 wind stresses. This wind stress field was used to 
determine how important preconditioning the model ocean with 
pre-El Nino wind stresses was in developing a realistic 
response. Figures 4 and 5 show the three month average wind 
stress fields used during spin-up of the early runs. Figures 
6-10 show the two year average wind stress fields used during 
the one year spin-up period of each of the different domain 
sizes. These figures show that, in general, these wind 
stresses decrease in strength from east to west and from high 
latitudes to low latitudes. This general pattern in the 
strength of the wind stress was also found in the analytic 
wind stresses used by Whitney (1992) in her final experimental 
runs. 

The last form of wind stresses used in this study were 
created by calculating the average wind stress for each month. 
These monthly values were interpolated to provide smoothed 
daily values. This was done by applying the monthly average 
wind stress to the 15th of the month. Values of the wind 
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stress between the 15ths of each month were found by linearly 
interpolating between the mid-month values. Since the data 
for December 1990 was not available, I applied the January 
1991 monthly average for the first 15 days of this month. The 
same technique was used to obtain values for the last 16 days 
of December 1992 . These time varying wind stresses based on 
the monthly averages were used to develop a baseline model run 
to determine the importance of daily wind stress variations. 

The 1991-1992 time period was chosen because of the well 
observed El Nino from the end of 1991 to May 1992, and the 
resurgence of the El Nino that took place at the end of 1992. 

C. DATA ANALYSIS TOOLS 

The analysis of the model data was greatly simplified with 
the use of the VIS5D graphics visualization software developed 
at the University of Wisconsin - Madison (Hibbard et al. 
1993) . Running this software on a Silicon Graphics 
workstation attached to a SGI 4D/380 VGX system, with 256 MB 
of RAM and 8 processors and graphics capable of displaying 10 6 
shaded polygons/sec, allowed for easy and rapid animation of 
all the model output fields using contours, vectors, and color 
shading. Animations of the model fields were used to 
determine when the model had reached a steady state during the 
spin-up stage. The determination of when the model had 
reached a steady state was done by animating the model's 
output fields and observing when the temperature and currents 
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had stopped changing with time. This software 
to verify the realism of the model predictions, 
invaluable tool in the analysis of the data, 
figures shown in this thesis were generated 
graphics tools. 



was also used 
VIS5D was an 
However, the 
with simpler 



19 



TABLE 1. MODEL PARAMETERS 



Parameter 


Symbol 


Value 


Units 


East-west grid spacing 


dx 


50/75/100 


km 


North south grid spacing 


dy 


25/var . 


km 


Time step 


dt 


600,1800 


s 


Heat flux coefficient 


Y 


6.9 x 10‘ 4 


Wm' 2o C' 1 


Horizontal eddy viscosity 


A 

m 


5.0 x 10 7 


cn^s" 1 


Vertical eddy viscosity 


k m 


1.0 


cm 2 s“ 1 


Horizontal eddy conductivity 


\ 


1.0 x 10 7 


cm 2 s" 1 


Vertical eddy conductivity 


k h 


1.0 


cm 2 s“ 1 


Thermal expansion coefficient 


a 


3.0 x 10- 4 


cm 2 s’ 1 


Reference density 


P 


1.0276 


gm cm* 3 | 
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Amin 


Am ax 


/ \ 


^ Q ^ +Tx f 




— <J>niaj 






f Ay 
















Level 


Depth (m) 








1 


10 








2 


20 








3 


30 








4 


50 








5 


80 








6 


110 








7 


150 








8 


200 








9 


300 








10 


500 








11 


800 








12 


1200 








13 


1700 


/ 




14 


2400 






15 


3500 


/ 









Figure 1. Illustration of the model domain coordinate 
system. 
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Figure 3. Temperature profile used to initialize the 
model . 
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Figure 5. Three month average wind stresses ( dynes /cm^ ) 
used to spin-up the model for the medium domain runs. 
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Figure 7 . Two year average wind stresses (dynes/cm^ ) usee 
to spin-up the model for the medium domain runs. 
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Figure 8. Two year average wind stresses 
(dynes/cm 2 ) used to spin-up the model for 
the large domain runs. 
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Figure 9. Two year average wind stresses 
(dynes/cm 2 ) used to spin-up the model for 
the large2 domain runs. 
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Figure 10. Two year average wind stresses (dynes/cnr) used 
to spin-up the model for the variable domain run. 
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III. MODELING RESULTS 



A. DEVELOPMENTAL MODEL RUNS 

Several developmental runs were performed to determine the 
dependence of the model’s response to variable wind stress 
forcing, basin size, grid size, and time step. Table 2 gives 
a listing of the different developmental runs performed. The 
major results from each of these runs are discussed in the 
following subsections in order to provide a basis for 
understanding the results from the experimental runs presented 
in the next section of this chapter. For most of the model 
runs, the surface (10 meter) temperature and currents, and a 
vertical cross section of the zonal velocity field along the 
equator are the only fields shown. These fields represent the 
major features of the model extremely well. 

1. Developmental run SI 

This model run used a three-month average wind stress 
field based on twice daily NOGAPS winds for December 1991 - 
February 1992 (Figure 4). The 1991-92 El Nino was near its 
mature phase during this period (Rasmusson and Carpenter 1982, 
Kousky 1992a) . The model was spun-up for 360 days until all 
the model fields appeared to have reached a steady state. 
Figure 11 shows the steady state currents and temperature 
structure in the upper (10 meters) level of the model. The 
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TABLE 2. DEVELOPMENTAL MODEL RUNS 
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strong meridional currents at the western boundary may be the 
result of the small basin size and the closed boundaries that 
prevent a realistic poleward outflow in the northwest and 
southwest corners of the domain. A cross section of the 
model's zonal velocity field along the equator (Figure 12) 
shows the vertical structure of the equatorial undercurrent 
(EUC). The model's EUC has a strong core in the eastern 
portion of the domain with the strongest flow having a maximum 
speed exceeding 20 cm/s at a depth between 50 and 80 meters. 
This is much less than the observed speed of the undercurrent 
(about 100-180 cm/s, Pond and Pickard 1983). Other models 
have also had problems in correctly representing the 
undercurrent speed maximum (e.g.. Miller et al. 1992). 

2 . Developmental Run S2 

The initial state for this model run was the final 
state of run SI. The model was forced with the daily varying 
wind stresses for December 1991 through February 1992 in an 
attempt to develop an El Nino-like response. The model’s 
initial response to the varying stresses was a slow 
strengthening of the cold tongue at the equator. Figure 13 
shows the temperature and currents in the upper level of the 
model at the time of maximum extent of the cold tongue. 
Figure 14 is a cross section at the equator of the model's 
zonal velocity field, at the same date as Figure 13. A 
comparison of Figures 11 and 13, and of Figures 12 and 14, 
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shows that the EUC slope has decreased but that the EUC still 
sholes to the west. Notice also that the EUC core has 
weakened slightly. Figures 15 and 16 show the same fields as 
Figures 13 and 14 at a time 24 days later when the cold tongue 
has retreated and an eastward surface current has developed 
between about 157 °W and 170 °W. The initial lengthening of the 
cold tongue during the first 40 days of the S2 run may 
indicate that some type of preconditioning in preparation for 
the retreat of the cold tongue is taking place. This tendency 
for the cold tongue to lengthen and narrow slightly can be 
seen in all of the model runs prior to a retreat of the cold 
tongue. Figures 17 and 18 from the end of the S2 run show a 
weak cold tongue at the end of the 91 day run. 

3. Developmental Runs S3 through S5 

The purpose of these runs was to investigate the 
sensitivity of the model to changes in the time step of 
integration. The time step was increased from an initial 
value of 600 seconds in increments of 1200 seconds until the 
model became unstable. At a time step of 3600 seconds the 
model was still stable and showed no signs of instability. At 
a time step of 4800 seconds the model was numerically 
unstable. No effort was made to determine the exact point at 
which the model would become unstable. 

The wind stress forcing for these runs was chosen so that 
the results could be compared to model run SI. Figures 19-24 
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show the temperature and currents for the model's upper level 
for each run, paired with cross sections along the equator of 
the model's zonal velocity field. Comparisons of Figures 19, 
21 and 23 with Figure 11 show that increasing the time step 
led to a general decrease in the 10 meter temperature, but a 
qualitatively similar temperature pattern. The 10 meter 
currents show little change with increasing time step. The 
cross sections of zonal currents (Figures 20, 22, and 24 
compared with Figure 12) show a slight decrease in the current 
speeds at all depths, however, this decrease is not enough to 
indicate that the time step tested should not be used. These 
results indicated that time steps ranging from 600 to 3 600 
seconds would be suitable for use in the later model runs. 

4. Developmental Runs S6 through S8 

Time steps of 1800, 2400, and 3600 seconds were used 
in model runs S6, S7, and S8 in order to ensure that the 
results obtained in model runs S3 through S5 would also apply 
to runs with variable wind stresses. These model runs are the 
same as run S2 but with larger time steps. Figures 25-30 show 
the temperature and currents for the model's upper level on 2 
February 1992 for runs S6, S7, and S8 paired with the 
corresponding cross sections along the equator of the model's 
zonal velocity field. Comparisons of Figures 25 and 27 with 
Figure 15 show that increasing the time step cooled the 10 
meter temperatures and strengthened the 10 meter currents. 
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Figures 26 and 28 show little difference in the zonal currents 
with depth compared to Figure 16. Figure 29 shows a small 
westward current in the midst of the eastward flow in the 
western portion of the domain. Figure 30 shows that this 
surface current is the manifestation of a large region of 
westward flow beneath the surface. 

Run S8 with a time step of 3600 seconds (Figures 23 
and 24) was stable. However, the region of westward currents 
that developed in the western portion of the domain (Figure 
24) was not observed in runs S6 or S7. This is probably an 
indication that under variable wind forcing this time step 
might be to close to the point of model instability. Even 
though Run S7 with a time step of 2400 seconds, (Figures 21 
and 22) gave results that compared reasonably well with runs 
S2 and S6 (time steps of 600 and 1800 seconds respectively) , 
a time step of 1800 seconds was chosen for subsequent runs to 
ensure that the model would remain stable under variable 
forcing. 

5. Developmental Runs Ml and M2 

Runs Ml and M2 were performed in order to determine 
the effects of basin size on the model's results. For these 
runs, the grid spacing used in the previous runs was 
maintained, but the number of grid points in both x and y was 
doubled, which quadrupled the area of the model domain. For 
run Ml the December 1991 through February 1992 three-month 
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average wind stress field (Figure 5) was used to spin the 
model up for 360 days. Due to the large region of cyclonic 
wind stress near the southwest boundaries of the domain, 
(Figure 5) the model was unstable for integration time steps 
greater than 600 seconds. Model run M2 used the final state 
of model run Ml as a initial state. In model run M2, the 
model became unstable when it was exposed to the time varying 
wind stresses no matter what time step was used. No attempts 
were made to correct the instability in this model run by 
relocating the domain, since the two-year average wind stress 
data was available, and it did not have this large cyclonic 
region. 

6. Developmental Run XI 

This model run was completed after experimental model 
runs M3 and M4 (described in Chapter III, Section B) had shown 
that use of the two-year average wind stress (Figure 7) 
corrected the instability problems of model runs Ml and M2. 
Based on the results of model runs M3 and M4 , another effort 
to increase the model's domain was made. Increasing the 
number of grid points in the model was not feasible due to 
computational run time limits. Therefore, the longitudinal 
grid spacing, ax, was doubled to 100 km, which increased the 
domain size to 120°x40°. This domain was large enough to 
study the remote effects of wind anomalies in the western 
Pacific on the eastern Pacific. The 1991-92 NOGAPS average 
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wind stress field was used to spin the model up for 390 days, 
at which time the model fields had reached a steady state. An 
error in the wind stress field used for the forcing of this 
run and run X2 was discovered near the end of model run X2 . 
This error caused the results in the extreme eastern portion 
of the domain to be incorrect, and it also increased the 2ax 
scale noise near the equator. However, the results for the 
rest of the domain were good and gave the information needed 
to determine the best position of the model domain for the 
final experimental run. 

Figure 31 shows the final temperature and current 
pattern in the upper level of model for this run. Figure 32 
shows a cross section of zonal velocity along the equator. 
There is a small area of weak eastward surface flow in the far 
western portion of the domain. This current is probably due 
to the extremely weak average winds at the equator in this 
region (Figure 8). Increasing the domain size also increased 
the strength of the equatorial undercurrent. 

Also of interest is the development of waves in the 
temperature field, which are most noticeable on the flanks of 
the cold tongue in the eastern Pacific. Figure 33 shows an 
example of the temperature waves along the middle of the cold 
tongue. The waves in this model run tend to be stronger along 
the northern flank of the cold tongue. This asymmetry is 
consistent with the observations of Legeckis (1977) and Harvey 
and Patzert (1976). The model's waves propagate to the west 
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at a speed of approximately 18-24 cm/sec, which is about half 
of the speed observed by Legeckis (1977). This difference may 
be due in part to the model's weak westward currents which are 
about half the speed of the observed currents. Semtner and 
Holland (1979) observed these waves in the steady state of 
their model but noted that the waves were too symmetric about 
the equator and were too strong in the southern hemisphere. 
This could be due to the constant wind stress field they used 
in contrast to the space and time varying field used in this 
model run. 

7 . Developmental Run X2 

Run X2 was started from the final state of run XI. The 
model was forced with the daily varying wind stress for the 
entire two year period. The equatorial cold tongue of the 
model retreated to the east three times during the two-year 
period. The first retreat of the cold tongue occurred between 
April 1991 and July 1991. This event was extremely small and 
was probably due to a few strong westerly wind bursts in the 
western Pacific. The second occurrence was in early November 
1991 and lasted for a period of 5 months. The third retreat 
was much weaker and occurred at the end of November 1992 and 
was still in progress at the end of the run. Figures 34 and 
36, respectively, show the temperature and currents associated 
with the greatest deviations from the initial state of the 
model for each of these events. Figures 35 and 37, 
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respectively, show the zonal currents in cross section along 
the equator for the same time periods. Figures 38 and 39 show 
these same fields at the point between the second and third 
retreats when the temperature and current structure had 
returned to almost the same State that was present at the 
beginning of the run. The second event showed a maximum 
retreat of the cold tongue of approximately 30° of longitude 
in late March 1992. In the third retreat, the cold tongue had 
moved only about 11° to the east by the end of the run. The 
wind stress fields for 1993 were not available to determine if 
the third event would have reached the same strength as the 
second event or if it would have remained weaker throughout 
its duration. Figures 35 and 37 show that not only did a 
eastward surface current develop during the retreats, but a 
westward equatorial undercurrent developed at about 100 meters 
in the western Pacific. This pattern appears to have split 
the equatorial region into two separate zonal circulations, 
one in the western Pacific and one in the central and eastern 
Pacific, instead of one large circulation across the entire 
basin (Figures 35 and 37). The two zonal circulations met at 
the surface in the vicinity of the dateline, where the surface 
horizontal flow was strongly convergent. 

B. EXPERIMENTAL MODEL RUNS 

Based on the results obtained in the developmental runs, 
three major experimental runs were performed. The three 
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domain sizes for these experiments were the ones labeled 
small, medium, and large2 in Figure 2. Table 3 lists the 
experimental runs in pairs, with the first run being the spin- 
up phase, and the second run being the run in which variable 
wind forcing was used. This section provides information and 
observations on all of the runs. However, the X3 and X4 runs, 
and comparisons with the X5 run, are the main focus of this 
section. The X3 and X4 runs gave the best simulation of the 
mean temperature and velocity fields and showed the most 
realistic responses to the variable wind forcing. 

1 . Experimental run S9 

The model for this run was spun-up from rest for 360 
days using the two year average wind stress forcing. The 
location of the model domain (Figure 2) was chosen so that 
some westerly wind bursts in the western Pacific would be 
contained in the model domain. Figure 40 shows the final 
temperatures and currents in the model’s upper level. Figure 
41 shows a cross section along the equator of the zonal 
velocity to show the equatorial undercurrent structure. The 
general structure of the fields is similar to that obtained by 
Semtner and Holland (1979), although this model shows less 
symmetry across the equator, probably due to the use of 
spatially and temporally varying wind stress fields. A 
comparison of this run (Figures 40 and 41) with run SI 
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(Figures 11 and 12) shows a colder and larger cold tongue and 
a more elongated EUC. 

2. Experimental run S10 

Starting from the final state of model run S9, the 
model was forced with the daily varying wind stresses for all 
of 1991-92. The model showed several realistic responses to 
these winds. The cold tongue retreated to the east on three 
different occasions. The first model retreat was a weak event 
that lasted from May 1991 through July 1991. During this 
retreat, the model showed a mild surface warming, probably in 
response to westerly wind stresses in the western and central 
Pacific early in 1991. Animations of the model’s temperature 
and velocity fields showed some indications of a Kelvin wave 
associated with this retreat of the cold tongue. Another 
study that used daily varying NOGAPS winds to force a 2 1/2 
layer ocean model has shown the passage of a Kelvin wave 
through the equatorial Pacific during the middle of 1991 
(Phoebus and Kindle 1993). In this study the model Kelvin 
wave had a very weak temperature and height anomaly signature 
compared to the Kelvin wave associated with the 1991-92 El 
Nino. According to the Climate Diagnostic Bulletin (Kousky 
1991), a weak warm event and negative values of the SOI were 
observed during this time frame. The real ocean and 
atmosphere didn't recover from this warm event as quickly as 
did the model. In the observations the warm event in mid 1991 
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blended directly into the large warm event at the end of 1991 
with only a small cooling and weakening of the SOI in between 
the two. 

The last two S10 model run retreats coincided with the 
observed retreats at the end of 1991 (Kousky 1992a) and again 
at the end of 1992 (Kousky 1992b) . Figure 42 shows the 
temperature and currents at the point of maximum retreat for 
the second retreat that started near the end of 1991. A cross 
section along the equator of the zonal velocity component 
(Figure 43) shows the corresponding equatorial undercurrent 
structure. In Figures 42 and 43 the reversal of the SEC 
associated with the cold tongue retreat has already 
dissipated. The EUC has reestablished itself across the 
entire domain but is tilted down from west to east, a reversal 
of its normal tilt. During this period, observations show the 
1991-92 El Nino became firmly established (Kousky 1992a) . 
However, in the model, the cold tongue tended to move rapidly 
back to its original position instead of remaining in a 
retreated state through May 1992, as was observed. Thus the 
westerly wind bursts that extended into this model's domain 
were capable of initiating an El Nino-like response in the 
model, but they were not able to sustain this pattern long 
enough. This suggests that westerly wind bursts to the west 
of the model domain (i.e. in the western Pacific) may be 
important to the development and maintenance of an El Nino. 
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Figures 44 and 45 show the S10 model run temperatures 
and currents at the end of the 1992. These figures represent 
the third cold tongue retreat that was developing in December 
1992. 

Before each of the three periods of retreat, the model 
SEC strengthens and the cold tongue tended to narrow and 
lengthen slightly. This is consistent with Wyrtki's (1975) 
hypothesis that the ocean must be preconditioned before an El 
Nino by the piling up of more warm water in the western 
Pacific. 

In the model's retreat of the cold tongue at the end 
of 1991, a strong reversal of the SEC occurred in the 
equatorial western Pacific. This strong eastward current had 
dissipated by the point of maximum retreat of the cold tongue, 
so it is not visible in Figures 42 and 43. Animations of the 
model's output fields showed that this eastward flow developed 
initially from a southward extension of the NECC. Below the 
surface the NECC also extended south and merged with the 
northern flank of the EUC. Then the NECC continued to extend 
equatorward as the EUC core sholed. This region of eastward 
flow then developed into a full reversal of the SEC in the 
equatorial western Pacific. Several aspects of this eastward 
equatorial flow can be seen in Figure 46, which shows the 
structure of the currents in a meridional cross section in the 
western portion of the domain on 13 November 1991, at the 
beginning of this merging of the currents. Figure 47 shows 



45 



the same field one day later when the currents have almost 
completely merged at the equator. This merging of the 
currents only occurred during the large retreat of the cold 
tongue at the end of 1991. Several strong reversals of the 
SEC in the western Pacific were observed during 1991-1992, 
including a reversal during June of 1992 that dissipated 
before the one associated with the 1992 El Nino (Kousky 
1992c) . 

Tropical instability waves like those observed in 
model runs XI and X2 were observed in this model run and in S9 
although the waves were much weaker. The model waves tended 
to weaken during the retreats of the cold tongue as seen in 
observations of these waves (Halpern 1988) . 

3. Experimental Run M3 

Run M3 was done to better investigate the effects of 
westerly wind bursts in the western Pacific and to determine 
the effects of the small basin size on the results obtained in 
runs S9 and S10. Run M3 is a repeat of run S9 but with the 
next larger basin size and a larger time step (cf . Table 3) . 
The longer time step was used to minimize the increase in run 
time arising from the larger domain. The developmental runs 
showed that although this time step increase might change the 
results quantitatively, it should not change the results 
qualitatively. The medium size domain model was spun-up from 
rest for 390 days to ensure that it had reached a steady 



46 



state. Figure 48 shows the temperature and currents of the 
model in its steady state. Figure 49 shows the zonal currents 
in cross section along the equator. Although the EUC is very 
weak and diffuse in the western portion of the domain, it is 
stronger and has a very well defined core over most of the 
domain. The model's EUC also developed a more realistic tilt 
in this run as compared to model run S9 (Figure 41), with the 
core of the undercurrent rising from west to east. The model 
also developed a narrower more well defined cold tongue than 
in model run S9 (Figure 40). 

4 . Experimental Run M4 

The final state of model run M3 was used as an initial 
state for model run M4 . This run is the larger domain and 
longer time step version of run S10. Due to the larger domain 
which extended further into the western Pacific, run M4 had 
more westerly wind bursts in the wind stresses than did run 
S10. In addition, the direct effects of the boundary 
conditions on the interior of the model basin were smaller. 
The M4 run showed the same cold tongue retreats seen in run 
S10. However, the retreat at the end of 1991, lasted about 2 
months longer, ending in March of 1992. By July 1992, the 
model had returned to a state resembling the steady state of 
run M3. 

Figures 50 and 52 show the temperature and currents at 
the point of greatest retreat for the second and third 
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retreats. Figures 51 and 53 show the zonal currents in a 
longitudinal cross section along the equator for the same 
periods. The model's upper level temperature tended to 
increase throughout this model run, indicating that perhaps a 
longer spin-up time was needed in model run M3 . Although the 
cold tongue does not appear to have retreated much in Figure 
52, the beginnings of an eastward surface current can be seen 
that would probably have led to a cold tongue retreat if the 
model run had been continued. 

The tropical instability waves seen in runs XI and X2 
were evident in this model run. The M4 run waves also had a 
strong seasonal cycle with the waves being strongest during 
August and September of 1991 and 1992 when the cold tongue was 
well developed. The waves could be seen throughout the model 
run; however, they were weaker during the strong cold tongue 
retreats at the ends of 1991 and 1992. 

5 . Experimental Run X3 

Based on the results of run X2 , the large model domain 
was altered slightly to include more area in the western 
Pacific. This put the western boundary at 135°E and the 
eastern boundary at 90°W. This alteration was done because an 
examination of the NOGAPS winds for 1991 and 1992 showed that 
a shift in the model domain further to the west (see Figure 2) 
would include almost every westerly wind burst that occurred 
in the Pacific during the period of interest. This westward 
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shift of the model domain also placed a region of generally 
weak wind stresses near the western boundary (Figure 9). This 
meant that the wind forcing in the western part of the domain 
was generally too weak to directly force strong currents. 
This tended to reduce the magnitude of the problematic western 
boundary currents and seemed to improve the model results for 
the entire western basin. After making these changes, model 
run X3 was spun-up from rest for 360 days, using the two year 
average wind stress field. At this time the model circulation 
had achieved a steady state. 

Figure 54 shows the temperature and currents in the 
upper level once model run X3 reached a steady state. Figure 
55 shows the zonal currents in a cross section along the 
equator. Figures 54 and 55 show an equatorial surface current 
in the western Pacific. The relaxation of the wind stress 
near the equator in the western portion of the domain yielded 
a negative wind stress curl which forced an eastward 
equatorial current. This current was much stronger than in 
the XI model run, probably due to the larger portion of the X3 
model domain in the western Pacific area of weak wind 
stresses. If a climatological value of the wind stress had 
been used, this eastward current would probably have developed 
more to the north at the latitude of the NECC. 

Figure 56 shows the model's mixed layer depth which is 
about 10 meters over most of the equatorial western Pacific, 
with the mixed layer depth getting deeper to the north and the 
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south. As discussed in Chapter II, the mixed layer depth is 
not allowed to go above the top model level (10 meters) . Thus 
the mixed layer depth in Figure 56 may be an overestimate of 
the actual depth of turbulent mixing. Figure 57 is a cross 
section along the equator of the model's temperature field. 
This Figure shows that the model's thermocline is between 80 
and 110 meters in the east and between 30 and 50 meters in the 
western portion of the domain. 

6. Experimental Run X4 

In run X4 , the final state of run X3 was used as a 
starting point from which the model was forced with the 1991 - 
1992 NOGAPS daily wind stresses. During the X4 run, the model 
ocean was exposed to strong westerly and easterly wind bursts 
in the western portion of the domain, and the eastward 
equatorial surface current seen in the western Pacific in run 
X3 (Figures 54 and 55) weakened and eventually disappeared 
(not shown) . Figure 58 shows the temperature and currents in 
the upper level of the model on 29 August 1991, after the 
eastward current on the equator had completely disappeared. 
Figure 59 shows the zonal currents in a cross section along 
the equator at the same time. 

The three retreats of the equatorial cold tongue seen 
in runs S10 and M4 were also seen in this model run although 
their structures were slightly different. The first and 
weakest retreat occurred before the eastward equatorial 
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surface current in the western Pacific had completely 
disappeared. Figures 60 and 61 show the upper level 
temperature and currents of the model at the point before this 
weak retreat began (26 March 1991) and at the time when this 
retreat was at a maximum (12 July 1991). 

Figures 62 shows the model's temperature and currents 
in the upper level on 16 October 1991. At this time, the 
equatorial cold tongue had reached its maximum extent. Figure 
63 shows the corresponding cross section of the zonal velocity 
component. This figure shows that, compared to run X3 (Figure 
55), the eastward equatorial surface current had diminished 
greatly, and the EUC had become better organized. Figure 64 
shows the model's mixed layer depth on 16 October 1991. 
Figure 65 shows a cross section of the temperature field along 
the equator for the same day. These figures provide a good 
basis for assessing the model's evolution during the first 
major retreat of the cold tongue. 

The second retreat of the model cold tongue began in 
the late October 1991. Figure 66 shows the temperature and 
currents on 11 December 1991 when a strong eastward equatorial 
surface current had formed in the western Pacific. Figure 67 
shows a cross section of the zonal currents along the equator 
for this same time. The source of this eastward current was 
similar to that in previous runs (eg. run S10). That is, the 
NECC in the western Pacific slowly extended southward toward 
the equator. The NECC then merged with the large reversal of 
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the SEC that developed in the model. In the spin-up run (X3) , 
and in this run the equatorial undercurrent did not extend 
into the western Pacific (Figure 67) . Thus, there was no 
strong eastward EUC in the western Pacific for the NECC to 
merge with, as in model run S10. If wind stresses from 1990 
were used to force the model, so that it was in a more 
realistic state at the start of the retreat of the cold 
tongue, a western Pacific EUC might have formed and merged 
with the NECC as in model run S10. Figures 68 and 69 show the 
zonal currents in a meridional cross section at 165 °E during 
the period when the NECC current was moving onto the equator. 

An unusual feature that developed in the model during 
the second cold tongue retreat was a tongue of warm water that 
extended eastward into the cold tongue (Figure 66) . This warm 
tongue was produced at the strong eastward equatorial surface 
current in the western and central Pacific rapidly advected 
warm water to the east along the equator. It is difficult to 
compare this feature to anything seen in the observations of 
the 1991-1992 El Nino, since the model does not actually 
represent the true surface temperature but a 10 meter 
temperature. I examined data from the TOGA buoys near the 
equator to see if they showed a comparable signal. However, 
the few buoys that did exist in the region where this feature 
was strongest in the model had several gaps in their data to 
provide a useful comparison. An examination of SST analyses 
from the U.S. Navy's optimal thermal interpolation scheme 
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(OTIS, Figure 70) shows that the SST tended to increase on the 
equator during the last 40 days of 1991, while the 
temperatures to the north and the south decreased slightly or 
stayed the same. Since the OTIS SST data is derived to a 
large extent from satellite data, surface heating and cooling 
may have masked the actual strength of the SST warming on the 
equator. The OTIS data is also spatially smoothed, and this 
would tend to mask the strength of this very narrow feature. 
This data does suggest that the model's eastward equatorial 
warm water intrusion may be realistic. However, the model may 
overdevelop the temperature gradients associated with the warm 
water intrusion, due to the extremely strong cold tongue that 
developed in the model . 

Figures 71 through 74 show the same model fields as 
Figures 62 through 65 but at the time of maximum retreat of 
the cold tongue on 13 March 1992. During the retreat of the 
cold tongue, the equatorial portion of the SEC had reversed in 
the western and central Pacific. The resulting strong 
eastward current advected warm water to the east. By mid- 
March 1992, this eastward current had begun to weaken, and the 
SEC was beginning to reestablish itself (Figures 71 and 72) . 

Although the restoration of the cold tongue in the 
model was mostly completed by 30 August 1992, the maximum 
extent of the cold tongue was not reached until 27 October 
1992 (not shown) . The cold tongue then began to retreat again. 
This third retreat was not nearly as strong as the second 
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event. At the end of 1992, the cold tongue had retreated 
about 10° of longitude, compared to a retreat of approximately 
20° at the end of 1991. Figure 75 shows the upper level 
temperature and currents on 31 December 1992 which was the end 
of run X4 . Figure 76 shows a cross section at the equator of 
the zonal velocity for the same time. The development of the 
cold tongue retreat seen in Figures 75 and 76 is comparable to 
that on 11 December 1991 (Figures 66 and 67), prior to the 
mature phase of the second cold tongue retreat. 

By using daily varying wind stresses instead of 
temporally averaged wind stresses, I was able to observe the 
model’s response to numerous typhoons that passed through the 
model domain. The expected response would be a swath of 
colder water in the region that the typhoon has just passed 
over (Ginis 1993). Observing the animated model output showed 
this effect in the wake of every typhoon that passed through 
the region. Although some of the typhoon cooling events were 
obscured by the large cold tongue in the model, several events 
were quite clear. 

Figures 77 through 82 show the wind stresses between 
the western boundary and 180 °E every two days during the 
period 13-23 November 1992. Typhoons Hunt and Gay can be seen 
intensifying and moving westward across this portion of the 
model domain. Figures 77-79 show that these typhoons are in 
a broad band of cyclonic flow in the northern hemisphere. In 
Figure 80, typhoon Hunt can be seen moving out of the area in 
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the upper right corner of the figure, and typhoon Gay 
continuing to move west northwestward across the domain. The 
most intense westerly winds on the equator in Figures 77-82 
are associated with typhoons Hunt and Gay and are part of the 
westerly wind bursts that forced the retreat of the cold 
tongue at the end of 1992. 

Figure 83 shows the X4 model run upper level 
temperatures and currents for 13 November 1992. Figures 84 
through 88 show the upper level temperature anomalies and 
currents for 15-23 November 1992, corresponding to the days 
represented in Figures 78-82. The temperature anomaly (TA) on 
a given day is calculated by subtracting the 10 meter 
temperature on 13 November 1992 from the 10 meter temperature 
on that day. For example the TA on 15 November 1992 is: 

•^lSAfov - -^lSAfov -^13 Nov 

Figures 84-88 show negative TA’s (i.e., cooling) in the wakes 
of typhoons Hunt and Gay. The TAs for 23 November 1992 
(Figure 88) show a large region of cooling in the wake of 
typhoon Gay of about 0.5° C with a maximum cooling of 2.5°C. 
The maximum cooling agrees well with observed temperature 
drops (e.g., 2-3°C, Ginis 1993), but should cover a larger 
region than is seen in the model. However, this model does 
not allow for factors such as reduced insolation due to clouds 
and the addition of fresh cool water in the form of 
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precipitation. The model's cooling is mainly due to mixing, 
upwelling, and advection. 

Figures 84-88 also show a region of cyclonic flow in 
typhoon Gay's wake. This cyclonic flow tends to be slightly 
west of the cooling region. This is consistent with cooling 
due to upwelling, with the maximum cooling lagging the maximum 
cyclonic flow by a few degrees of longitude and a day or two. 
Figure 89 shows a north-south cross section through this 
divergent region on 21 November 1992. This figure shows that 
the model produced strong upwelling in this region down to 500 
meters, with an extremely strong signature down to 150 meters. 

The model ' s mixed layer depth deepened with the 
storm's passage (not shown). The region of colder water in 
the wake of the typhoon slowly returned to normal but was 
still visible 10 days after the storm (not shown) . This 
strong response to the realistic wind forcing shows the 
ability of the model to handle real winds and the importance 
of realistic daily varying wind stresses in forcing the ocean. 

The tropical instability waves on the equator were 
quite distinct in model run X4 (not shown). The seasonal 
cycle to these waves developed as in runs X2 , S10, and M4 . 
The waves were observed throughout the duration of the model 
run, although their intensity and the number of waves varied 
with the season. The waves were not as strong in late 1991 
and early 1992 when the cold tongue was quite weak. The waves 
were stronger during July-September , when the NECC and SEC 
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were relatively strong and the cold tongue had extended well 
into the central Pacific. The waves propagated to the west at 
about 20 cm/s, which is about half the observed speed 
(Legeckis 1977). The model's weaker and slower propagation 
speed may be due to the weak SEC and NECC which had about half 
the observed speed (12 cm/s versus the observed 25 cm/s for 
the NECC, and 40 cm/s versus the observed 60 cm/s for the SEC 
Pickard and Emery 1990). This would result in a weaker shear 
environment in which to develop the waves and a slower 
westward motion of the waves. 

7 . Experimental Run X5 

Model run X5 was performed to provide a comparison to 
model run X4 that would highlight the role of the daily 
varying winds. Run X5 repeated X4 but used smoothed wind 
stresses instead of the daily fluctuating wind stresses. 
Chapter II, section B, page 17, describes the procedure used 
to obtain these smoothed wind stresses. 

This model run was started late, so the results are 
only available through the end of 1991. The model’s response 
to the smoothed winds is similar to that of model run X4 in 
general structure and response. However, several differences 
developed in run X5. As expected, these differences are 
largest at later dates in the model runs. The absence of 
strong typhoons in run X5 kept the temperatures warmer in 
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general, probably due to the lack of the intense mixing and 
upwelling associated with these wind events. 

Figures 90 and 91 show the upper level temperature and 
currents on 26 March 1991 and 12 July 1991. Comparing runs X4 
and X5 (Figures 90 and 60) shows that the basic current 
patterns are the same, but many of the small scale current 
variations are missing from the smoothed run. This is 
consistent with the lack of rapid variations in the wind 
stress. At the later date, a comparison of runs X4 and X5 
(Figures 91 and 61) shows that although no westerly winds 
existed in the smoothed fields for the summer of 1991, the 
easterly winds did relax enough for the model to develop the 
weak retreat of the cold tongue that had been observed in 
previous runs. 

Figure 92 shows temperature anomalies calculated by 
subtracting the fields of model run X4 from those of X5 for 1 
October 1991 (just before a retreat of the cold tongue). This 
figure shows that the western equatorial Pacific is generally 
warmer, and the eastern equatorial Pacific is generally cooler 
in model run X5 than in X4 . The warmer west Pacific and 
cooler east Pacific regions in run X5 may be due to a 
reduction in eastward warm water advection from the west 
Pacific resulting from a reduction in westerly wind bursts in 
run X5. The cooler eastern Pacific in run X5 also may be due 
to increased upwelling caused by the more constant easterly 
winds resulting from the smoothing of westerly wind bursts in 
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the run X5 wind stresses. Figure 93 shows a cross section 
along the equator of the run X5 temperature anomalies. This 
figure shows that the strong anomalies are limited to the 
upper 150 meters of the model domain. 

Figure 94 shows the upper level current anomalies for 
run X5 on 1 October 1991. The differences between the current 
patterns of the two models is small, but numerous small 
anomalies do exist. The mixed layer depth of model run X5 is 
on average 5-10 meters shallower than that of model run X4 , 
due to a decrease in the turbulent mixing (not shown) . The 
anomaly fields became quite complex for the period 
corresponding to the retreat of the cold tongue (not shown) . 
This is because a difference in position of the cold tongue of 
a degree of longitude is enough to generate extremely large 
temperature and current anomalies. 

Figure 95 shows the upper level temperatures and 
currents for 11 December 1991. Figure 96 shows a cross 
section along the equator of the zonal velocity field. These 
figures can be compared with Figures 66 and 67 which show the 
same fields on 11 December 1991 in model run X4. Comparisons 
of these figures show similar major features exist in both 
model runs. The SEC has reversed, and warm water has been 
advected to the east. The large scale changes in the model 
tend to dominate the figures, and so the small scale 
differences in the two models are not as apparent as they are 
when the cold tongue is well established. The speeds of the 
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major currents (EUC, SEC, and reversed SEC) are all faster in 
this model run than in X4 . This is probably due to the more 
constant wind forcing throughout model run X5. The cold 
tongue also retreated slightly farther to the east by the end 
of 1991 in run X5 than in X4 . This is again probably due to 
the weaker, but more constant, westerly wind forcing that 
existed near the equator when all of the westerly wind bursts 
were smoothed together. 

In general, although the results of model run X5 were 
similar to those of X4 , numerous differences developed due to 
the lack of the daily varying winds. The major current speeds 
were 10-20 cm/s higher when using smoothed wind stresses 
instead of daily varying wind stresses. Based on the results 
obtained by the end of 1991, the differences between the two 
model runs would probably continue growing through the rest of 
the X5 run, with large differences developing during the 
retreat of the cold tongue at the end of 1992. 

C. SUMMARY OF RESULTS 

The model's response to the variable wind forcing used in 
this study was excellent. The model remained stable and 
developed realistic responses to wind forcing events on time 
scales of a few days to years. The increases in basin size 
allowed the model to develop more realistic current and 
temperature patterns. This also reduced the problems due to 
large anomalous currents that tended to develop on the western 
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boundary in the small basin runs, and reduced problems with 
wave reflections off the model's closed boundaries. The daily 
varying wind field was extremely important in the development 
of the model's temperature and current patterns. Experimental 
run X5 shows that the absence of daily variations in wind 
stress may lead to large basin scale temperature and current 
anomalies . 

One area of concern in the early developmental work on 
this model (Whitney 1992) was the extremely noisy mixed layer 
depth field. The increase of the basin size has tended to 
reduce the development of pockets of extremely shallow or deep 
mixed layer depths. Some of the remaining small scale 
variations in the mixed layer depth may be explained by the 
mixed layer depth definition used in this model. The mixed 
layer depth is defined by the depth of turbulent mixing, which 
can have large variations over short horizontal scales due to 
changes in surface currents and wind forcing. 
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model run. Dashed (solid) contours indicate westward 
(eastward) flow. 



62 




currents (cm/s) on 9 January 1992, showing the maximum 
extent of the cold tongue. 
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currents (cm/s) on 2 February 1992, when the cold tongue 
has retreated to the east. 




zonal currents (cm/s) along the equator on 2 February 
1992. Dashed (solid) contours indicate westward (eastward) 
flow. 
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model run (29 February 1992). Dashed (solid) contours 
indicate westward (eastward) flow. 
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1800s. 




Figure 20. Depth-longitude cross section of model run S3 
zonal currents (cm/s) along the equator at the end of the 
model run with a At of 1800s. Dashed (solid) contours 
indicate westward (eastward) flow. 
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Figure 21. Model S4 upper level temperature (°C) and 
currents (cm/s) at the end of the model run with a At of 
2400s. 




Figure 22. Depth-longitude cross section of model run S4 
zonal currents (cm/s) along the equator at the end of the 
model run with a At of 2400s. Dashed (solid) contours 
indicate westward (eastward) flow. 
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currents (cm/s) at the end of the model run with a At of 
3600s . 




Figure 24. Depth-longitude cross section of model run S5 
zonal currents (cm/s) along the equator at the end of the 
model run with a At of 3600s. Dashed (solid) contours 
indicate westward (eastward) flow. 
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zonal currents (cm/s) along the equator on 2 February 
1992, with a At of 1800s. Dashed (solid) contours indicate 
westward (eastward) flow. 
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currents (cm/s) on 2 February 1992, with a At of 2400s. 




Figure 28. Depth-longitude cross section of model run S7 
zonal currents (cm/s) along the equator on 2 February 
1992, with a At of 2400s. Dashed (solid) contours indicate 
westward (eastward) flow. 
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westward (eastward) flow. 
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model run. Dashed (solid) contours indicate westward 
(eastward) flow. 
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south boundaries of the eastern portion of the cold tongue 
at the end of model run. 
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currents (cm/s) on 21 March 1992. During the second cold 
tongue retreat . 




Figure 35. Depth-longitude cross section of model run X2 
zonal currents (cm/s) along the equator on 21 March 1992. 
During the second cold tongue retreat. Dashed (solid) 
contours indicate westward (eastward) flow. 
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cold tongue retreat. 




1992. During the third cold tongue retreat. Dashed (solid) 
contours indicate westward (eastward) flow. 
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(eastward) flow. 
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model run. Dashed (solid) contours indicate westward 
(eastward) flow. 
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cold tongue retreat. 




1992. During the second cold tongue retreat. Dashed 
(solid) contours indicate westward (eastward) flow. 
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tongue retreat. 




zonal currents (cm/s) along the equator on 31 December 
1992. During the third cold tongue retreat. Dashed 
(solid) contours indicate westward (eastward) flow. 
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Figure 46. Meridional cross section of run S10 U velocity 
at 166°W on 13 November 1991, showing the merging of the 
NECC and EUC. Dashed (solid) contours indicate westward 
(eastward) flow. North is to the right. 




(solid) contours indicate westward (eastward) flow. North 
is to the right. 
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currents (cm/s) at the end of the model run. 




model run. Dashed (solid) contours indicate westward 
(eastward) flow. 
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During the second cold tongue retreat. Dashed (solid) 
contours indicate westward (eastward) flow. 



82 




currents (cm/s) on 31 December 1992. During the third cold 
' tongue retreat. 




Figure 53. Depth-longitude cross section of model run M4 
zonal currents (cm/s) along the equator on 31 December 
1992. During the third cold tongue retreat. Dashed (solid) 
contours indicate westward (eastward) flow. 
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Figure 54. Model run X3 upper level temperature (°C) and 
currents (cm/s) at the end of the model run. 




model run. Dashed (solid) contours indicate westward 
(eastward) flow. 
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end of the model run. 




(°C) along the equator at the end of the model run. 
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Figure 59. Depth-longitude cross section of model run X4 
zonal currents (cm/s) along the equator on 29 August 1991, 
once the equatorial eastward surface current has 
disappeared. Dashed (solid) contours indicate westward 
(eastward) flow. 
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reached it's maximum extent before the first retreat. 




currents (cm/s) on 12 July 1991. During the first cold 
tongue retreat. 
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zonal currents (cm/s) along the equator on 16 October 
1991. Time of greatest extent of the cold tongue. Dashed 
(solid) contours indicate westward (eastward) flow. 
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1991. Time of greatest extent of the cold tongue. 
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eastward surface current in the equatorial western 
Pacific . 




1991. Note the eastward surface current in the equatorial 
western Pacific. Dashed (solid) contours indicate westward 
(eastward) flow. 
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zonal currents (cm/s) at 165°E on 13 November 1991. Note 
the eastward surface current north of the equator. Dashed 
(solid) contours indicate westward (eastward) flow. 




current onto the equator. Dashed (solid) contours indicate 
westward (eastward) flow. 
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Figure 70. OTIS SST time series for the last 40 days of 
1991 at three stations bisecting the equatorial warm 
tongue seen in model run X4 (Figure 66). 
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tongue retreat. 




Figure 72. Depth-longitude cross section of model run X4 
zonal currents (cm/s) along the equator on 13 March 1991. 
During the second cold tongue retreat. Dashed (solid) 
contours indicate westward (eastward) flow. 
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Figure 73. Model run X4 mixed layer depth in meters on 13 
March 1992. During the second cold tongue retreat. 




Figure 74. Depth-longitude cross section of model run X4 
temperature field (°C) along the equator on 13 March 1992. 
During the second cold tongue retreat. 
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tongue retreat . 




zonal currents (cm/s) along the equator on 31 December 
1992. During the third cold tongue retreat. Dashed (solid) 
contours indicate westward (eastward) flow. 
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Figure 77. Wind stresses ( dynes /cm^ ) for 13 November 1992. 
Note weak typhoon Hunt at about 158°E, 8°N and the early 
stages of typhoon Gay located on the eastern edge of the 
figure at about 7°N. 
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(warmer) temperatures. Contour interval is .5°C. 
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November 1992. Dashed (solid) contours indicate cooler 
(warmer) temperatures. Contour interval is .5°C. 




(warmer) temperatures. Contour interval is .5°C. 
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November 1992. Dashed (solid) contours indicate cooler 
(warmer) temperatures. Contour interval is .5°C. 




(warmer) temperatures. Contour interval is .5°C. 
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field (cm/s) at 159°E on 2 1 November 92, in the region of 
upwelling near typhoon Gay. Vertical velocity exaggerated 
by factor of 86400. Longest vector equals vertical 
velocity of .2 m/day 
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(°C) on 1 October 1991. Positive anomalies indicate model 
run X5 was warmer than X4 . Dashed (solid) contours 
indicate cooler (warmer) temperatures. Contour interval is 
. 3 °C . 




contours indicate cooler (warmer) temperatures. Contour 
interval is .5°C. 
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October 1991. 
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currents (cm/s) on 11 December 1991. Showing the strong 
eastward surface current on the equator. Compare to run X4 
( Figure 66 ) . 




Figure 96. Depth-longitude cross section of model run X5 
zonal currents (cm/s) along the equator on 11 December 
1991. Compare to run X4 (Figure 67). Dashed (solid) 
contours indicate westward (eastward) flow. 
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IV. SUMMARY 



A. DISCUSSION 

In this study, a tropical ocean general circulation model 
with an imbedded mixed layer has been used to investigate the 
tropical Pacific ocean's response to daily varying winds and 
westerly wind bursts. 

The model was tested using several different basin sizes 
with the larger basin sizes giving the most realistic results. 
The smaller basin size runs did develop the same major 
features as the larger runs and were very useful for studying 
small scale or short time frame effects. The increases in 
basin size led to faster current speeds, with the EUC reaching 
a maximum speed in run X4 of greater than 70 cm/s. This is 
still slower than the observed speeds of approximately 150 
cm/s (Pond and Pickard 1983) but is a large improvement over 
the speeds (about 30 cm/s) seen in the small basin runs. 

The runs with the largest domain size (runs X3, X4 , and 
X5 ) and with two years of daily varying wind stress for 1991- 
1992 developed three retreats of the cold tongue. The first 
retreat was a weak event compared to the later two retreats. 
This weak retreat may have been associated with passage of a 
weak Kelvin wave across the model domain that was created by 
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strong westerly wind bursts in March through May of 1991 (cf. 
Phoebus and Kindle 1993). 

The second major retreat of the cold tongue corresponds 
well with the mature phase in December 1991 - May 1992 of the 
1991-92 El Nino. This retreat can again be tied to an 
increase in westerly wind bursts in the western Pacific during 
September - November of 1991. Westerly wind bursts continued 
to occur near the equator in the western Pacific into early 
1992. These westerly winds bursts provided the forcing needed 
to advect warm water to the east. 

The third retreat of the model ' s cold tongue came at the 
end of 1992 and coincided with the resurgence of the El Nino. 
This last retreat started later and developed slower than the 
one a year earlier. Wind stress data from 1993 was not 
available to see if this equatorial warming in the model at 
the end of 1992 would have developed as the observed El Nino 
did. 

The model also developed several other features that 
correspond well with observed phenomena in the Pacific ocean. 
The model developed tropical instability waves along the 
equator during the spin-up phase of the model runs. These 
waves were also present in the variable wind forced runs, but 
varied in intensity with the strength of the NECC and the 
SECC. The size and seasonality of these waves match well with 
similar observed waves. 
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The model’s response to individual westerly wind bursts 
was checked by looking at the model’s response to the strong 
typhoons that moved through the model's domain. The model 
developed regions of decreased upper level temperatures in the 
wake of the typhoons due to increased mixing and upwelling. 

The model quite clearly shows the importance of the 
westerly wind bursts that are associated with tropical 
cyclones that form near the equator. The cold tongue retreat 
at the end of 1991 was associated with several typhoons that 
moved through the western Pacific during October and November. 
In particular typhoons Yuri and Tia in November 1991 provided 
the major westerly wind burst forcing for the the model’s cold 
tongue retreat. In November 1992, typhoons Hunt and Gay 
produced extremely strong westerly winds on the equator that 
appear to have initiated the retreat of the model’s cold 
tongue at the end of 1992. 

Comparisons of model runs X4 and X5 showed that the daily 
variations in the wind stresses used for this study played an 
important role in developing both the small scale and 
transient features, and in developing the basin wide 
temperature patterns. The temperature differences between 
these two runs showed that the smoothed winds tended to cause 
more upwelling and advection of cold water to the west in the 
eastern half of the model domain, leading to colder eastern 
Pacific SST's. This is due to the lack of westerly wind 
bursts in the eastern Pacific that would tend to decrease the 
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upwelling or advect warm water to the east. The smoothed 
winds also tended to reduce the amount of warm water advected 
out of the western Pacific leading to higher SST's in this 
region. Again, this is caused by the lack of westerly wind 
bursts that would advect the water to the east. Only during 
the El Nino period do the westerly wind bursts become strong 
enough and persistant enough to show up in the average fields 
as westerly winds. 

B. CONCLUSIONS 

The model's response to the inclusion of the daily 
fluctuating winds from NOGAPS gave very encouraging results 
and allowed the model to develop a realistic El Nino-like 
response. Comparisons of the model runs forced with daily 
varying wind stresses and those forced with smoothed wind 
stresses showed that these differences in forcing produced 
significant short term, small scale differences in the upper 
ocean currents and temperatures. In addition, these forcing 
differences led, over several months of model integration, to 
intraseasonal , basin wide differences in the currents and 
temperatures. The runs with daily varying wind stresses also 
indicate that individual tropical cyclone forcing events play 
a major role in the model's development of cold tongue 
retreats and warming of the central and eastern Pacific. 
Thus, the model's simulation of the 1991-1992 El Nino was 
enhanced by the use of realistic daily varying wind stresses. 
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The model's response to the winds is even more impressive 
given the many simple forcing parameterizations used (e.g., 
the surface heat flux) . For example, the model was able to 
develop a retreat of the cold tongue without complex 
temperature and salinity forcing. This highlights the role of 
westerly wind bursts in the development of an El Nino-like 
response. 

The model's development of an El Nino-like response also 
shows the usefulness of this model in the study of this major 
climatological event. The NOGAPS wind forcing data could be 
obtained and applied to the model on a monthly basis to verify 
early indications of an El Nino. 

Mismatches between the model's upper level temperatures 
and the observed SST pattern in the ocean are probably the 
major flaws in the model's depiction of the tropical Pacific. 
The development in the model of a more realistic subtropical 
gyre circulation for the return of water to the eastern 
Pacific should make the long narrow cold tongue in the model 
more closely resemble the observations. 

The model did not properly resolve the deep currents and 
jets that have been observed in the tropics due to the lack of 
detailed resolution in the deep ocean. The model also did not 
develop realistic western boundary currents due to the lack of 
bottom topography and realistic coastlines. The proper 
modeling of these features would improve the character of the 
model results. However, their inclusion should not change the 
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major features observed in the model's response. This is in 
partdue to the efforts that were made in this study to reduce 
the effect either one of these features would have on the 
upper levels in the central Pacific which is the primary area 
of focus in this study. 

C . RECOMMENDATIONS 

Several options for improving the model results are being 
investigated, and additional efforts to improve the model 
domain are needed. The inclusion of realistic coastlines and 
bottom topography should improve the model's results. The 
north and south boundaries need to be moved further from the 
equator so the north and south subtropical gyres are properly 
represented. This increase in the domain should improve the 
model's upper level temperature field in comparison to 
observations . 

However, the model used for runs X3, X4 , and X5 can still 
be used to investigate many phenomena associated with El Nino. 
The inclusion of a spatially and temporally varying heating 
parameter would allow the determination of how variable 
heating modifies the model's simulation of El Nino. 
Comparison of the model's results with and without this 
heating would give an indication of the relative importance of 
the wind and heating in developing and maintaining an El Nino. 

Model runs using monthly climatological winds interpolated 
to daily values need to be done to determine baseline 
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variations of the model fields. A comparison of these fields 
with the results of model run X5 would highlight the 
importance of the differences in the average wind field in El 
Nino years versus non-El Nino years. This model run would 
also provide a more realistic spin-up run for use in later 
model runs. 

Several efforts were made to include biharmonic diffusion 
in the model to replace the current Laplacian diffusion. The 
model is unstable when run with both mixed layer physics and 
biharmonic diffusion. Further work in this area is needed to 
determine a method of using biharmonic friction with the mixed 

layer physics. 

* 

I was able to identify individual tropical cyclones that 
played an important part in initiating and maintaining the 
cold tongue retreat. However, no effort has been made to try 
to identify specific events that lead to the preconditioning 
of the cold tongue prior to these events. The identification 
of individual events that led to the stronger easterlies that 
created the proper conditions for the cold tongue retreat 
would be a significant step forward in understanding El Nino. 

The model output also needs to be thoroughly analyzed 
using the same type of statistical methods applied to the 
observational data by Cooper (1992) and Kent (1993). This 
work would provide a better way of comparing the model data to 
the observations than just using visual comparisons. This 
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would allow 
patterns and 



t more careful verification of the model's wave 
other temporal and spatial variations. 
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